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Abstract 



Using a nonperturbative approach we examine the large frequency asymp- 
totics of the two-point level density correlator in weakly disordered metallic 
grains. This allows us to study the behavior of the two- level structure factor 
close to the Heisenberg time. We find that the singularities (present for ran- 
dom matrix ensembles) are washed out in a grain with a finite conductance. 
The results are nonuniversal (they depend on the shape of the grain and on 
its conductance), though they suggest a generalization for any system with 
finite Heisenberg time. 

PACS numbers:71.30.-Fh, 05.60.-Fw, 72.15.Rn 
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A great variety of physical systems are known to exhibit quantum chaos. The common 
examples are atomic nuclei, Rydberg atoms in a strong magnetic field, electrons in disordered 
metals, etc Chaotic behavior manifests itself in the energy level statistics. It was a 
remarkable discovery of Wigner and Dyson, that these statistics in a particular system 
can be approximated by those of an ensemble of random matrices (RM). Here we consider 
deviations from the RM theory taking an ensemble of weakly disordered metallic grains with 
a finite conductance g as an example. The results seem to be extendible to general chaotic 
systems. 

There are two characteristic energy scales associated with a particular system: a classical 
one Ef. and a quantum one. The quantum energy scale is the mean level spacing A. In a 
chaotic billiard, for example, is set by the frequency of the shortest periodic orbit. Well 
developed chaotic behavior can take place only if E'c ^ A. 

In a disordered metallic grain the classical energy is the Thouless energy E^. = D/L^, 
where D is the diffusion constant, and L is the system size. For a weakly disordered grain 
the two scales are separated by the dimensionless conductance g = Ec/ A ^10. For 
frequencies u <^ Ec the behavior of the system becomes universal (independent of particular 
parameters of the system ). In this regime in the zeroth approximation the level statistics 
depend only on the symmetry of the system and are described by one of the RM ensembles: 
unitary, orthogonal or symplectic 0. 

One of the conventional statistical spectral characteristics is the two-point level density 
correlator 

Kiu, x) = {pie + u,H + x6H)p{e, H)) - A'^ (1) 

where H is the Hamiltonian of the system, 5H is a perturbation, x is the dimensionless 
perturbation strength and p(e, H + x5H) = Tr5 (e — H — x6H) is the x-dependent density of 
states at energy e. It is convenient to introduce the dimensionless frequency s = u/A and the 
dimensionless two-level correlator R{s,x) = A'^K{uj,x). Dyson determined R{s,x = 0) 
for RM. For example, R{s, o) in the unitary case equals to 

, ^, , 1 — cos(27rs) , ^ 

R{s,0) = S{s) ^ ' 2 

and is plotted in the insert in Fig. |T[ 

Perhaps the most striking signature of the Wigner-Dyson statistics is the rigidity of 
the energy spectrum Among the major consequences of this phenomenon are: a) the 
probability to find two levels separated by c<j ^ A vanishes as u ^ 0; h) the level number 
variance in an energy strip of width A^A is proportional to InA^ rather than A^; and c) 
oscillations in the correlator R{s, 0) in Eq. (Q) decay only algebraically. 

In the two level structure factor |^ S{t,x) = J^^dsexp{iTs)R{s,x) the reduced fiuc- 
tuations of the level number manifest themselves in the vanishing of S{t, 0) at r = 0, and 
the algebraic decay of the oscillations in R{s, 0) leads to the singularity in S{t, 0) at the 
Heisenberg time r = 27r. In the unitary case, e. g. 

5(r,0) =min{|r|/(27r),l}. 
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At r ^ 27r this Dyson result was obtained by Berry for a generic chaotic system by use of 
a semiclassical approximation. To the best of our knowledge nobody succeeded in analyzing 
the behavior of S{t, 0) around t = 2n using this formalism. 

Wigner-Dyson statistics become exact in the limit g = Ec/A oo. We consider correc- 
tions to these statistics for finite g. One of the better understood systems in this respect is 
a weakly disordered metallic grain. For frequencies much smaller than Ec the statistics are 
close to universal ones, the corrections being small as {s/gY @1- At s ^ 1 the monotonic 
part of R{s, x) can be obtained perturbatively p 

i?p(s,a:)=3?5:^^-^^— ^, (3) 

where are eigenvalues (in units of A) of the diffusion equation in the grain, a = 2 for the 
unitary ensemble and a = 1 for the orthogonal and symplectic ensembles . At this point 
we can define 

E, = eiA/vr^, g = ejir^, (4) 

where ei is the smallest nonzero eigenvalue. Perturbation theory allows one to determine 
S{t, 0) at small times r <^ 1. Since the oscillatory part of R{s,x) is non-analytic in 1/s it 
can not be obtained perturbatively. 

In this Letter we obtain the leading s 3> 1 asymptotics of R{s, x) retaining the oscillatory 
terms and monitor how the singularity in S{t, 0) at the Heisenberg time is modified by the 
finite conductance g. We make use of the nonterturbative approach |TT| that is valid for 
arbitrary relation between s and g. The oscillatory part Rosc{s,x) = R{s,x) — Rp{s,x) for 
the unitary (u), orthogonal (o) and symplectic (s) cases equals to 



RUs,x) = '-^§^P{s,x), (5a) 



R:Js,x) = -'-^^P\s,x\ (5b) 



ns / N cos(7rs) ^, , cos(27rs) n \ 

RUs,x) = -^-^P{s,x) - -^^P\s,x), (5c) 

where j/^ = + x^, and P{s, x) is the spectral determinant of the diffusion operator 

1 

(6) 



Pis,x)= n 



Note that Eq. (|]) expresses Rp{s, x) through the Green function of this operator. Thus, 
regardless of the spectrum e^, Rp{s,x) and Rosc{s,x) are related: 

It follows from Eq. (^) that -Rosc("S, x) together with P(s, x) decays exponentially ai s g. 
As a result, the singularity in S{t, 0) at the Heisenberg time is washed out: S{t, 0) becomes 
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analytic around r = 27r. The scale of smoothening of the singularity is 1/Ec { see Fig. p. At 
1 -C s -C 5^ the sum of Eq. (H) and Eq. (H) gives the leading high frequency asymptotics of 
the universal result, fors^g it coincides with the perturbative result Rp{s,x) of Ref. 0. 

In a closed d-dimensional cubic sample (diffusion equation with Dirichlet boundary con- 
ditions) = gir'^n'^, where n = [rii, . . . ,nd) and rii are non-negative integers. For s ^ g 
and (i < 4 we obtain the asymptotics P(s, 0) — > exp{— 7r(s/7r(7)°'/^/[r((i/2)(isin(7r(i/4)]}. At 
1 ^ s <^ (7 we obtain 



sm^fvrs) sin^fTTS) 



, , oil J. \/i a I oiii \/i 'J J ^ — ^ X /o\ 

R[s,0) = — \ / 2-2^2' (8) 

This result was shown in Ref. to be valid even for s < 1. Thus, it is natural to assume 
that for the unitary ensemble the sum of Eq. ( |Sa]) and Eq. (^ gives the correct g —>■ oo 
asymptotics at arbitrary frequency. This is related to the absence of higher order corrections 
to the leading term of the perturbation theory S{t, 0) oc r in the unitary case. 

Now we sketch the derivation of our results. Consider a quantum particle moving in 
a random potential V{r). The perturbation acting on the system is a change in the po- 
tential 6V{r). Both V{f) and 5V^(f^ are taken to be white noise random potentials with 
variances {V{r}V{r')) = 6{f- r')/2'nvT and {5V {r)5V {r')) = x^A6{f- P)/(47rz/), At < 1, 
{V{r)6V{r')) = 0, where () denotes ensemble averaging and z/ is the density of states per 
unit volume. The dimensionless perturbation strength x'^ is assumed to be of order unity. 



We use the supersymmetric nonlinear a- model introduced by Efetov [|r^ , and follow his 
notations everywhere. One can show that for the system under consideration the a-model 
expression for K{uj, x) is given by 



1 

Kicu,x) = / VQeM-Fj 



(9a) 



j=o 



Til' 

~8 



J df^Ti <D{VQf + 2iujAQ + iJAkQ - ^(AQ)4 . (9b) 



The 8x8 supermatrix Q{f) obeys the nonlinear constraint Q{rY = 1 and takes on its values 
on a certain symmetric space H = G/K, where G and K are groups ||T^. For example, for 



the unitary ensemble H = SU(1, l/2)/SU(l/l) SU(1/1) [|13|. The integration measure 



for Q in the functional integral Eq. (Pa|) is the invariant measure on H 



The hierarchy of blocks of supermatrices is as follows: advanced-retarded (A-R) 
blocks, fermion-boson (F-B) blocks, and blocks corresponding to time-reversal. A = 
diag{l, 1,1,1,— 1,-1,— 1,-1} is the matrix breaking the symmetry in the advanced- 
retarded (A-R ) space, k = diag{l, 1,-1,— 1,1,1,— 1,-1} is the symmetry breaking matrix 
in the Fermion-Boson ( F-B ) space. 

The large frequency asymptotics of K{uj, x) can be obtained from Eq. (pa]) by use of the 
stationary phase method. The conventional perturbation theory corresponds to integrating 



over the small fluctuations of Q around A [11 




Q = Ail+tP)il-tP)-\P= " ^ , (10) 
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where the matrix P describes these small fluctuations. 

Q = A is not the only stationary point on H. This fact to the best of our knowledge was 
not appreciated in the literature. The existence of other stationary points makes the basis 
for our main results. 

It is possible to parameterize fluctuations around a point Qq in the form Q = Qq{1 + 
iPo) (1— zPo)~^- Expanding the Free Energy Fj in Eq. (|9b| ) in Pq we would obtain the station- 
arity condition dFj/dPo = 0. This route however is inconvenient because the parametriza- 
tion of Pq will depend on Qq. Instead we perform a global coordinate transformation on 
H that maps Qq to A, Qq — >■ T(^^QqTq = A. We note that the matrices A and —Ak be- 
long to H, and the corresponding terms in Eq. (^) can be viewed as symmetry breaking 
sources. This transformation changes the sources, but allows us to keep the parametrization 
of Eq. ( p!OD and preserves the invariant measure. Introducing the notation Q\ = TQ~^Aro 
and Q\k = TQ^AkTo we rewrite Eq. (j^) as 

2 

Kiu,x) = '^^IvqI dfj rfr-'STr(QAA:g(rl)Str(gAfcg(r-^))exp(-F(QA)), 
F{Qa) = — j rfr-STr | D{VQf + 2iujQ^Q - -^{QkQY \ . (11) 

The stationarity condition dF{Q;C) / dP\p^Q = implies that all the elements of Q\ in 
the AR and RA blocks should vanish (this can be seen from Eq. (p!OD). 

Here we discuss in detail only the calculation for the unitary ensemble. The calculation 
for the other cases proceeds analogously, and we just point out the important differences 
from the unitary case. 

Consider now the unitary case. The 4x4 supermatrices B and B in Eq. ( PUD are given 

by 

-as \ 

a 

i6* ' 
iai ih ) 

where a, h and their conjugates are ordinary variables, and ai, a2 and their conjugates are 
grassmann variables. 

The only matrix besides A that satisfles the stationarity condition is Qa = A = —kA. 
In this case Q\k = —A. All other matrices from H contain nonzero elements in the AR and 
RA blocks. Both stationary points contribute substantially to K{uj,x). 

Consider the contribution of Qa = A to K{u!,x) flrst. We substitute Qa = —kA and 
Q^k = —A into Eq. (|TT]) and expand both the Free Energy F{Q/^) and the pre-exponent to 
the second order in B and B. Expanding B{r) in the eigenf unctions of the diffusion operator 
0^(r) as B{r) = J^n 4>fj.{^Bfj_, and introducing Sfj, = is + + x^ we arrive at the following 
expression for the dimensionless density-density correlator: 

Rl,is,x) = ^f VB,{Y.{K\' + \b,\' - ai>i,, + a2>2,^])' 

X exp(-27r{-zs + + - - <,^2,,)]})- (12) 



B 
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We have to keep the perturbation strength finite to avoid the divergence of the inte- 
gral over caused by the presence of the infinitesimal imaginary part in uj. For the non 
parametric case we should take the — limit after the integral in Eq. ( |T^ ) is evaluated. 

Since the Free Energy in Eq. (|l^) contains no Grassmann variables in the zero mode they 
have to come from the pre-exponent. Therefore out of the whole square of the sum in the 
pre-exponent only the terms containing all four zero mode Grassmann variables contribute. 
In these terms the prefactor does not contain any variables from non-zero modes. Thus, 
the evaluation of the Gaussian integrals over non-zero modes yields the superdeterminant 
of the quadratic form in the exponent. Supersymmetry around A is broken by s, therefore 
this superdeterminant differs from unity and is given by P{s,x) of Eq. (^). The correctly 
ordered integration measure for Grassmann variables is Y{^daifj,dal^dal^da2fj.- Evaluating 
the integral we arrive at Eq. (|5a|). 

In quasi-lD for closed boundary conditions and x = the spectral determinant 0) 
can be evaluated exactly, and from Eq. (^) we obtain 

^ sinhM^/f +sin2 



For Qa = A the same procedure as used above leads to Eq. (|^), which coincides with the 
result of Ref. g. 

The behavior of ^(r, 0) at r = and r = 27r is associated respectively with i?p(s,0) 
(Eq. (H)) and R^^^{s,Gi) (Eq. (|5aD ). In other words the singularity at the Heisenberg time 
is determined by the contribution to R{s, 0) from A. It is clear that the cusp in S{t, 0) at 
r = 27r will be rounded off because i?Qg^(s,0) decays exponentially at large s. The scale of 
the smoothening is of order 1/ g. 

The Fourier transform of Eq. (p!3[) ( see Fig. ) is 

Qu.o^,,n^ ^ (-irexp(-7r^n^^|t|) |t| 
^ 7r"'(yfnsmh(7ra) 47r 

Even though S'"^(27r + 1, appears to be a function of \t\, it is regular at t = 0. 

We can also estimate S'"(27r,0)^ in any dimension. It is proportional to 1/g of Eq. (^ 
and is given by 

5'^(27r,0); ' ■ 

Consider now T-invariant systems. For the orthogonal ensemble there are still only two 
stationary points on H: A and A = —kK. To determine the contribution of the A-point 
we use the formula Eq. (|lT]) with Qa = A and QkK = —A and Efetov's parametrization 
for the perturbation theory [llj. The calculations are analogous to those for the unitary 




ensemble and lead to Eq. (pBI). The contribution of Qa = A gives Eq. (^. At r = 27r the 
third derivative of S{t, 0) for the orthogonal ensemble has a jump. This singularity also 
disappears at finite g. 
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In the symplectic case there are three types of stationary points which correspond to 
singularities in the structure factor S{t,0) at r = 0,7r, 27r §]. The r = 27i singularity 
corresponds to Q\ = A = —kA, and its contribution to R{s,x), given by the second term 
in Eq. is exactly the same as i?°g^(s,x). The stationary point Q\ = A corresponds to 
the r = singularity in S{t,0) and leads to Eq. The r = tt singularity corresponds to 
a degenerate manifold of matrices Qa on H Q\ = diag(r^, U2, — r^, — 12), QkA = —kQA, 
where IL2 is a 2 x 2 unit matrix, = nixTx + ^y'Ty, = 1 and r^^y are Pauli matrices in 
the time-reversal block. The calculation proceeds as before and leads to the first term in 
Eq. (^. In quasi-lD we can obtain the leading contribution to the structure factor S{t, 0) 
around r = tt 

f°° — 4:8110^ (q\t\z)dz , / X ^/ , X / X 



The result is plotted in Fig. 0. In all dimensions the logarithmic divergence in the zero mode 
result is now cut off by finite g, and S^{n, 0) oc In^f. 

In conclusion we mention several points about our results. 1) Equation (|) describes the 
deviation of the level statistics of a weakly disordered chaotic grain from the universal ones. 
This deviation is controlled by the diffusion operator. This operator is purely classical. It 
seems plausible that the nonuniversal part of spectral statistics of any chaotic system can 
be expressed through a spectral determinant of some classical system-specific operator. If 
so, the relation Eq. (0) should be universally correct! 

2) The formalism used here should be applicable even to the systems weakly coupled 
to the outside world (say through tunnel contacts). As long as the level broadening F 
(u; = ^uj + zF) is smaller than Ax^ the integration over the zero mode variables in Eq. (p!2D is 
convergent. The integral over the other modes is always convergent provided F < E^. Thus, 
the presence of a perturbation can effectively "close" a weakly coupled system. Under these 
conditions Eq. (^ remains vahd after the substitution cos(27rs) — » exp(— 27rF/A) cos(27rs) 
and a;^ — > — F/A. 

3) The classification of physical systems into the three universality classes (unitary, or- 
thogonal and symplectic) is, of course, an oversimplification. In practice there is always a 
time scale which determines the crossover from one ensemble to another. For example if a 
system is subjected to a magnetic field for very short times it will still effectively remain 
orthogonal. On the other hand, the long time behavior will be unitary. The characteristic 
time is set by the strength of the magnetic field. 

For a disordered metallic grain in a magnetic field this characteristic time is given by 
I'h/D. For frequencies larger than D/ljj the system effectively becomes orthogonal. This 
implies that even if we neglect the spatially nonuniform fluctuations of the Q-matrix the 
cusp in S{t, 0) at r = 2tt will be washed out on the scale of Alj^/D ( although there will still 
remain a jump in the third derivative of S{t, 0) ). For the system to behave as unitary for 
frequencies of order Ec the magnetic length Ih has to be shorter than the size of the system. 
Spin-orbit interaction that causes the orthogonal-to-symplectic crossover can be considered 
in a similar way. 

4) The rounding off of the singularity in S'(27r, 0) is also present in the random matrix 
model with preferred basis [jl4[ Note that our results differ from those in Ref. |15 
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substantially. This means that finite g is not equivalent to finite temperature for the corre- 
sponding Calogero- Sutherland model |16|. 

We are grateful to D. E. Khmel'nitskii, B. D. Simons and N. Taniguchi for numerous 
discussions throughout the course of this work. 
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FIGURES 




FIG. 1. Structure factor in quasi- ID case for unitary symmetry (solid line) and the universal 
structure factor (dashed line). Insert: the two level correlator as a function of level separation. 
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FIG. 2. The structure factor for the symplectic case in quasi-lD (sohd hne) and the universal 
result (dashed line). Insert: the universal stucture factor. 
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